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ABSTRACT 

We present the results of CO ( J = 3 — 2) and HCO+ ( J = 4 — 3) mapping observations toward a 
nearby embedded cluster, Serpens South, using the ASTE 10 m telescope. Our CO ( J = 3 — 2) map 
reveals that many outflows are crowded in the dense cluster-forming clump that can be recognized as a 
HCO+ clump with a size of ~ 0.2 pc and mass of ~ 80 Mq. The clump contains several subfragments 
with sizes of ~ 0.05 pc. By comparing the CO ( J = 3 — 2) map with the 1.1 mm dust continuum 
image taken by AzTEC on ASTE, we find that the spatial extents of the outflow lobes are sometimes 
anti-correlated with the distribution of the dense gas and some of the outflow lobes apparently collide 
with the dense gas. The total outflow mass, momentum, and energy are estimated at 0.6 Mq, 8 Mq 
km s~^, and 64 Mq \vo? s~^, respectively. The energy injection rate due to the outflows is comparable 
to the turbulence dissipation rate in the clump, implying that the protostellar outflows can maintain 
the supersonic turbulence in this region. The total outflow energy seems only about 10 percent the 
clump gravitational energy. We conclude that the current outflow activity is not enough to destroy 
the whole cluster-forming clump, and therefore star formation is likely to continue for several or many 
local dynamical times. 

Subject headings: ISM: clouds — ISM: jets and outflows — stars: formation — submillimeter — 
turbulence 



1. INTRODUCTION 

Most sta rs form in clusters ()Lada fc Ladal 120031 : 

lAUen et "aD 120061 ). Therefore, understanding the for- 
mation process of star clusters is a key step towards 
a full understanding of how stars form. Recent obser- 
vations have revealed that star clusters form in turbu- 
lent, m agnetized, parsec- scale dense clumps of molecular 
clouds (jRidge et al.ll2003l ) . These clumps contain masses 
of 10^ — lO'^ Mq, fragmenting into an assem bly of cores 
that coll apse to prod uce stars (e.g. .Davis et al.l 119991: 



SandcU fc iSIei [200TI: ?: lAndre et all 120071: iKirk et al 
2007i: ,Walsh"et al.ll2007l: ISaito et al. l l2008l:lMaru'ta' et al 



20101 ). In cluster- forming clumps, stellar feedback 
such as protostellar outflows, stellar winds, and radi- 
ation rapidly start to shape the surroundings. Be- 
cause of the short separations between forming stars 
and cores, these feedback mechanisms are expected to 
contr o l subsequent sta r forrn a tion (e.g., iNorman fc Silkl 



19801 iKru mholz et al.' [2001 iMcKee fc Ostrike"? 



Matzneri 120 07: Nakani urafc Lil I2007t IPeters et all 
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2010; 



Wang et al.ll2010l) . However, the roles of the stellar feed- 
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back on cluster formation remain poorly understood ob- 
servationally. 

According to recent theoretical studies, two main sce- 
narios of cluster formation are now under debate. For 
these two scenarios, the roles of the stellar feedback on 
cluster formation are different. One is the rapid, dynam- 
ical formation model, in which the large-scale turbulent 
flow mainly controls star formation in cluster-forming 
clumps, and star formation is envi s ioned to complete in a 
short time (e.g., lElmegreenl 120071 : iHartmann fc BurkertI 
[200l . In this case, cluster formation is expected to be 
terminated by stellar feedback due to initial star burst 
within a couple of turbulence-crossing times, so that star 
formation efficiencies (SFEs) can stay low at a observed 
level in nearby cluster-forming regions (SEE ~ 0.1). To 
promote rapid gravitational collapse and thus rapid star 
formation, the magnetic fields should be weak dynam- 
ically. Another scenario is the quasi- virial equilibrium 
model, in which cluster formation is considered to con- 
tinue at least for sev e ral dynamical, or free -fall, times 
(e.g., iTan et al.l [200l iLi fc Nakammal [20061 ). Because 
supersonic turbulence dissipates rapidly in a turbulence- 
crossing time, additional turbulent motions should be 
injected to maintain the supersonic turbulence. Pro- 
tostellar outflow feedback is proposed to pl ay a sig- 
nificant role in the t urbulence regenerat i on ([Matzneii 
[20071: INakamura fcH ^Mft ICarroll et al.1 12010D . The 
moderately-strong magnetic fields are also important to 
impede the rapid global gravitational collapse. In this 
case, the global inflow and outflow are expected to coex- 
ist, interacting with themselves. As a result, the cluster- 
forming clumps as a whole can keep quasi-equilibrium 
states for a relatively long time. 

In an effort to discriminate between the two scenar- 
ios and constrain the theory of cluster formation, we 
have made a start on a systematic observational study 
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of a nearby infrared dark cloud (IRDC), Serpens South, 
with a focus on the stellar feedback. Recent observa- 
tions strongly suggest that IRDCs are i n the very early 
phase of star cluster f ormati on ( Rathb orne et al.l 120061 : 
iPerreto k FviieMOO^^ iButler fc Tanll2009| ) and presum- 
ably still retain the primordial structure of the cluster's 
natal cloud. In particular, in the early phase of evolu- 
tion, the protostellar outflows are likely to play a more 
dominant role in cloud dynamics, compared to the stellar 
winds and radiation, because they are more powerful in 
the early protostellar evolution (Bontcmps et al. 1996). 
Therefore, IRDCs are expected to be suitable to explore 
the roles of protostellar outflows on cluster formation. 

Serpens So uth is a nearby ernbedde d cluster, recently 
discovered by iGutermuth et al.l (|2008[ ) using the Spitzer 
Space telescope. The accurate distance to Serpens South 
remains unknown. Most of the previous studies adopt 
the same distance as the Serpens Cloud Core, located 
about 3° north, because the YSOs associated with Ser- 
pens South have the same LSR velocities as the Serpens 
Cloud Core, whose distance is usually adopted as 260 ± 
37 pc (jStraizvs et al.llT9960 . However, larger distances 
of 415±25 pc, 380 pc, and 3601^3 pc have been recently 
claimed for the Serpens Cloud Core based on a VLBA 
parallax of EC95, a young AeBe s tar embedded in Ser- 
pens Cloud Core (IDzib et al.l lMol) SEDs of YSOs bas e 
on the near infrared observations ("Gorlova e t al.ll2010[ ). 
and X-ray luminosity func tions ( j Winston e t al.||2010D , re- 
spectively. According to IGutermuth et all (j'2008). the 
Serpens South cloud appears to be seen in absorption 
against PAH emission from the W40 HII region. There- 
fore, the distance to W40 is likely to be an upper limit to 
the distance to Serpens South, although even the former 
has not yet been deter mined to any satisfactory preci- 
sion (300 pc - 900 pc; IRodnev fc Reipurth 2008, '). In 
this paper, we assume 260 pc as the distance to Serpens 
South. 

From the Spitzer IRAC observations, IGutermuth et al.l 
(|2008[ ) revealed that the cluster contains about 60 
Class I protostars and 40 Class II sources. Recently, 
IBontemps et al.l (|2010D have discovered 7 Class can- 
didates in t his region through the Herschel Gould Belt 
Survey (see lAndre et al.l I2010D . In the central part of 
the cluster, the fraction of Class I protostars relative 
to all the YSOs (Class I/II) reaches about 80 % at the 
high surface density of 430 pc~^. These observations in- 
dicate that Serpens South is undoubtfully in the very 
early stage of cluster formation, i.e.. Serpens South is 
a "protocluster" that are forming a number of stars. 
IGutermuth et all ((2008) also estimated the median pro- 
jected distance between nearest neighbor YSOs to be 
about 3700 AU, signiflcantly shorter than the typical 
length of protostellar outflow lobes, implying that the 
outflows can potentially influence star formation in this 
region. 

Recent dust continuum observations have revealed that 
the cluster is located at the constricted region in a long 
filamentary cloud that is recognized as a dark cloud 
in mid-infrared (^ lOAtm) wave length (jGutermuth et all 
l200l. I20ITI: I Andre et all I20I0I ). Several sub-filaments 
also appear to converg e toward the main filament (e.g. 
Mvers 2009; Bontcm ps et al.l l2010t ). Very recently, 
Sugitani et al.l (|201l[ ) carried out near-infrared (JHKs) 



polarization observations toward Serpens South using 
the imaging polarimeter SIRPOL which is a polarimetry 
mode of the near-infrared camera SIRIUS mounted on 
the IRSF 1.4 m telescope at the South Africa Astronom- 
ical Observatory. They found that the global magnetic 
field is spatially well-ordered and almost perpendicular 
to the main filament. Several sub-filaments appear to be 
parallel to the global magnetic field. These observational 
facts imply that moderately-strong magnetic field played 
a role in the course of the formation of the Serpens South 
filament. The moderately-strong magnetic field appears 
consistent with the quasi-virial equilibrium model where 
the magnetic fields signifi cantly impede the glob al gravi- 
tational collapse (see, e.g.. lNakamura &:Till201lD . In this 
paper, we investigate the outflow activity in the Serpens 
South cloud on the basis of the CO (J = 3 — 2) and 
HCO+ ( J = 4 - 3) observations using the ASTE 10 m 
telescope. 

The rest of the paper is organized as follows. First, we 
describe the details of our CO ( J = 3 — 2) and HCO"*" 
( J = 4 — 3) observations toward Serpens South in Sec- 
tion [21 We present in Section [3] our CO ( J = 3 — 2) 
and HCO"'" (J = 4 — 3) maps toward Serpens South, 
which reveal a number of powerful coUimatcd outflows 
and clumpy structure in the central part of the protoclus- 
ter, respectively. We also compare our CO ( J = 3 — 2) 
map with the three-color Spitzer IRAC image (3.6 fim in 
blue, 4.5 fj,ia in green, 8.0 /im in red) to help identify the 
molecular outflow lobes. In the three-color image, the 
knots created by shocks due to the outflows are found as 
exten ded green objects w ith excess in 4.5 /xm emission 
(e.g.. lTakami et al.(l20ldl ). We derive the global proper- 
ties of the outflows identified from our CO ( J = 3 — 2) 
map. Then, in Section |4l we briefly discuss how the ob- 
served outflows influence the clump dynamics. Finally, 
we summarize our main conclusion in Section [5l 

2. OBSERVATIONS 

The CO (J = 3-2; 345.796GHz) and HCO+ (J = 4-3; 
356.734GHz) mapping ob servations were car ried out with 
the ASTE 10 m telescope (jEzawa et al.ll2004l ) toward Ser- 
pens South during the period of August 12 — 14 2010, 
using the on-the-fly (OTF) mapping technique (Sawada 
et al. 2008). Figure [T] presents the observation boxes 
that are overlaid on the col umn density map d erived from 
the Herschel observations ([Andre et al.ll2010l ) , indicating 
that our observation area covers most of the dense flla- 
mentary structure associated with Serpens South. Our 
main observation box (Box 1 in Fig. [T]) covers a 13' x 13' 
region centered at the dense part of the protocluster, 
(RA, Dec) = (18h30m3.4s, -2°2'4.6"). We make the 
image, by combining scans along the two axes that run 
at right angles to one another using the PLAIT algo- 
rithm developed by Emerson & Graeve (1988), in order 
to minimize the "scanning effect" . The observation re- 
gion was extended to a larger 20' x 20' box (Box 2 in Fig. 
[11). We also observed a 13' x 13' box (Box 3 in Fig. [J) 
centered on (RA, Dec) = (18h29ml6s, -1°41'35"). The 
extended observation areas (Boxes 2 and 3) have higher 
noises and suffer from the scanning effects due to the 
OTF observations because only one scan was done. The 
beam size in HPBW of the ASTE telescope is 22", which 
corresponds to 0.028 pc at a distance of 260 pc. The 
main-beam efficiency, 77, was 0.57 at 345 GHz. We used 
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a 345 GHz SIS heterodyne receiver, which had the typical 
system noise temperature of about 200 K in DSB mode 
at the observed elevation. The temperature scale was de- 
termined by the chopper-wheel method, which provides 
us with the antenna temperature, T^, corrected for the 
atmospheric attenuation. As a backend, we used four 
sets of 1024 channel autocorrelators, providing us with a 
frequency resolution of 125 KHz that corresponds to 0.11 
km s""'^ at 345 GHz. After subtracting linear baselines, 
the data were convolved with a Gaussian-tapered Bessel 
function (Magnum et al. 2007) and were resampled onto 
a 5" grid. The resultant effective FWHM resolution is 
24". The typical rms noise levels are 0.11 K and 0.51 K 
in with a velocity resolution of 0.5 km s^^ for Box 1 
and Box 2 and 3, respectively. 

As for the Spitzer IRAC data, archival data in three 
IRAC bands (3.6, 4.5, and 8.0 ^m) were obtained for 
Serpens South. All the data had been reduced with 
the basic calibration data (BCD) pipelines developed by 
Infrared Processing and Analysis Center (IPAC), and 
mosaicked with the Mosaicking and Point Source Ex- 
traction (MOPEX) software developed by Spitzer Sci- 
ence Center. The mean FWHMs of the point response 
functions (PRFs) are 1".66, 1".72, 1".98, for the three 
bands, respectively (See Spitzer Space Telescope Ob- 
server's Manual Version 8.0). 

3. RESULTS 
3.1. Global Cloud Properties 

Since both the CO (J = 3 - 2) and HCO+ ( J = 4 - 3) 
transitions have higher critical densities (approximately, 
10* and 10^ cm^'^, respectively) and higher transition 
energies (~ 33 and 43 K, respectively) than those of 
the lower transition lines such as CO {J — 1 — 0), it 
allows us to examine higher density and/or higher tem- 
perature gas. In FigureHl we present the CO ( J = 3 — 2) 
velocity integrated intensity map in the ranges from 
^^LSR = —2 km s~^ to -f 15 km s~^, which shows that the 
CO ( J = 3 — 2) emission tends to be much stronger in 
the southern box with about 20' x 20' area (Boxes 1 and 
2), while a relatively weak emission can be seen in the 
northern box with about 13' x 13' area (Box 3). In fact, 
the Serpens South protocluster resides in Box 1, and the 
strong C0(J = 3 — 2) emission originates mainly from 
the powerful molecular outflows of the cluster member 
YSOs, as described in the next subsection. We note that 
the weak CO ( J = 3 — 2) emission was detected in the 
entire observed area, indicating that the molecular gas is 
widely distributed in much larger area. 

We also present in Figures [3^, [S^, and [St, the CO 
(J = 3 - 2) peak intensity, HCO+ (J = 4 - 3) peak 
intensity, and 1.1 mm continuum maps, respectively, to- 
ward Box 2. The 1.1 mm continuu m map was taken by 
the AzTEC ca mera (iWilson et all [20 08) on the ASTE 
10 m telescope (jGutermuth et al.ll2011l). The positions 
of Cla ss I and II sources identified by iGutermuth et al.l 
(j2008f ) from the Spitzer IRAC observations are plotted 
in Figures [3^ through [31;. The 1.1 mm map indicates 
that the dense gas is distributed along a long filament 
running from north to south. This filament (hereafter, 
the main filament) has several sub-filaments that appear 
to converge toward the main filament. Such structures 
appear to be common in many star-forming region, as 



pointed out bv lMversI (|2009D . 

The CO (J = 3 — 2) peak intensity map in Fig. [2^ 
is not dominated by emission at the systemic velocity of 
the cloud, ~ 7.5 km s~^. Such a feature is prominent 
in the velocity channel maps presented in Figs. [4| and 
[5l where the emission is weak in the range from 7 
to ^ 10 km s^^. For comparison, we also present in 
Figures M and [7| the CO (J = 3 - 2) and HCO+ (J = 
4 — 3) spectra, respectively. The CO ( J = 3 — 2) spectra 
often show double-peak line profiles, whereas the HCO"*" 
( J = 4 — 3) spectra often show single-peak line profiles. 
By comparing the CO (J = 3 - 2) and HCO+ (J = 
4 — 3) spectra, it is clear that the CO ( J = 3 — 2) line 
profiles appear to be caused by self-absorption in the 
range from ~ 7 to '-^ 10 km s~^, rather than multiple 
velocity components. 

Strong CO (J = 3 - 2) and HCO+ ( J = 4 - 3) emis- 
sion come from almost the same part where the 1.1 mm 
emission is strongest and the Serpens South protoclus- 
ter is located. In fact, the CO (J = 3 - 2) and HCO+ 
( J = 4 — 3) emission peaks coincide reasonably well with 
the 1.1 mm peak. The HCO"*" ( J = 4 — 3) emission is 
distributed in an elongated clump with the axis ratio of 
about 1.5-2. The HCO+ (J = 4-3) integrated intensity 
map presented in Fig. [7| suggests that the clump con- 
tains several subfragments. The clumpy structure of the 
HCO+ clump is more evident in the HCO+ ( J = 4 - 3) 
velocity channel maps presented in Fig. [8l The typical 
size of the subfragments is about 0.05 pc. The positions 
of two subfragments (sources 1 and 2 indicated in Fig. 
[Tib) are indicated by the diamonds in Fig. [Sj As shown 
in the next subsection, these two subfragments contain 
the possible driving sources of the outflows. The HCO+ 
( J = 4 — 3) emission takes its maximum at the position 
of source 2 presented in Fig. [Tl] that coincides with the 
position of the 1.1 mm peak. The positions of these two 
subfragments coincide reasonably well with those of the 
3 mm compact sources recently identified by the NMA 
observations (Fukuda et al. in prep.). 

The CO line profiles have prominent blueshifted 
and/or redshifted wings near the region where the CO 
emission is strongest, as can be seen in Fig. [6l The 
line wings sometimes extend over the velocity range of 
— 10 ~ 30 km s~^. These prominent wings indicate the 
presence of powerful molecular outflows in this region. 
The CO ( J = 3 — 2) peak intensity map indicates that 
there are several local peaks and ridges near the central 
dense part of the protocluster. As shown in Section [3721 
these peaks trace the outflow lobes reasonably well. Sim- 
ilar high- velocity wings can also been seen in the HCO"^ 
( J = 4 — 3) line profiles (see Fig. [7]) that are are some- 
times extended in the range from about 3 km s~^ to 15 
km around the subfragments located near the 1.1 
mm peak. 

In the southern part of Box 2, the redshifted compo- 
nents of the identified outflows tend to be extremely weak 
(see also the line profiles presented in Fig. [61d). It is 
unclear what causes the extremely-asymmetric line pro- 
files in the southern area. One possibility is the exis- 
tence of a foreground cloud in the velocity of around 10 
km s~^ in the southern area. The foreground cloud, if 
cold, preferentially absorbs the redshifted outflow com- 
ponents. Another possibility is the existence of the global 
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infall toward the southern part of the cloud. If the cloud 
envelope is infalling, the foreground infall gas having ve- 
locities similar to the redshifted outflow components can 
absorb the redshifted components. In the next subsec- 
tion, we identify the molecular outflow lobes from CO 
( J = 3 — 2) data and discuss several features of the out- 
flow lobes in detail. 

3.2. Molecular Outflows 

Here, using the CO (J = 3 — 2) data cubes, we at- 
tempt to identify high-velocity components that origi- 
nate from molecular outflows. First, we scrutinize the 
velocity channel maps (Figs. |4]and[5]) and the position- 
velocity diagrams (see e.g., Fig. [T7|) to find localized 
blueshifted or redshifted emission. As can be seen in 
Figs. |3]and[5l the ^^CO high velocity structure is very 
complex and crowded in the southern box with about 
20' X 20' area (Box 2). Thus, identifying outflows from 
the ^^CO emission alone is extremely difficult in this 
crowded area. Therefore, we compare the CO image 
with three-color Spitzer IRAC images (3.6 ^m in blue, 
4.5 /zm in green, 8.0 /im in red) of this region, which 
are presented in Figure 1101 The blow-up of the cen- 
tral area is s hown in Figure 1111 According to previous 
studi es (e.g., iNoriega-Crespo et al.l \20Q^ iTakami et al.l 
l2010f ). extended objects that have strong emission in 
IRAC channel 2 (4.5 /zm) compared to IRAC channel 1 
(3.6 /im) and stand out as extended green objects in the 
three-color images are likely to be the objects shocked 
by the protostellar outflows. Here, we refer to such ob- 
jects as infrared H-H objects. To identify the infrared 
H-H objects, we first made an image of the remaining 
emission in the Spitzer IRAC band 2 (4.5 iim) image af- 
ter the subtraction of the emission in band 1 (3.6 /im). 
This image processing in practice allows us to discrimi- 
nate between shocks and stars, displayin g them with pos- 
itive and negative values , respectively ( Zhang fc Wangj 
I2009t ITakami et al.|[2010l ). We then visually inspected 
the image to search for mid-infrared outflows. We se- 
lected spatially-extended sources that have significant ex- 
cess at 4.5 fim as mid-infrared outflows. The positions of 
the identifled infrared H-H objects are indicated by the 
crosses in Figs. [TUb. fTTb. andfT^. Some of the infrared 
H-H objects have bow-shapes. For example, the infrared 
H-H objects labeled as Kl and K2 in Fig. [T2k have arc- 
like features, both of which appear to move toward the 
south west direction (see Figs. [13k and [l2b) . From our 
CO (J = 3 - 2) and Spitzer IRAC data, we identifled 15 
blueshifted and 10 redshifted components. The positions 
of the identified lobes are indicated in Fig. [SI Table [1] 
gives a brief summary of the identified outflow lobes. We 
also show in Fig. [THIthe blueshifted and redshifted high- 
velocity components identifled from HCO"*" (J = 4 — 3). 
In the following, we describe some of the main features 
of the high-velocity components we identifled from the 
CO (J = 3 - 2) data. 

The relatively strong blueshifted components Bl and 
B4 are well localized in the channel maps and appear 
at the maps within the wide range of 0.75 km s~^ to 
5.75 km s~^. These lobes have extremely high-velocity 
blueshifted components. In particular, the Bl compo- 
nent has moderately-strong emission even in the veloc- 
ity range of about —10 km s~^ (see Fig. [6]). There are 
also many infrared H-H objects associated with these two 



components. The Bl and B4 appear to move away from 
the central dense part toward north and south, respec- 
tively. In fact, the most distant infrared H-H object as- 
sociated with B4 has a bow shape, consistent with the 
interpretation that it is moving away from the dense part 
(see Fig. [TTb). As mentioned below, the B4 component 
might originate from a compact 3 mm source located at 
the 1.1 mm peak, recently identifled by the Nobeyama 
Millimeter Array (NMA) observations (Fukuda et al. in 
prep.). If this is the case, the corresponding redshifted 
components are likely to be R2 and R3. The position of 
the 1.1 mm peak is indicated by a diamond (source 2 in 
Fig. Illb ) in the Spitzer image. It appears that no clear 
Spitzer sources are associated with the 1.1 mm peak. The 
3 mm source discovered by the NMA observations may 
be a deeply-embedded, extremely-young protostar such 
as a Class source. In fact, the HCO+ ( J = 4 — 3) line 
proflle toward source 2 shows blue-shifted part that are 
stronger than the red-shifted part (see Fig. [H|) . Such a 
blue-skewed proflle of an optically-thick line may be in- 
dicative of infalling gas, suggesting that the 3 mm source 
located at the 1.1 mm peak is an infalling object having 
a powerful outflow. 

There are two faint blueshifted components labeled 
with B2 and B3. They are vaguely seen in the channel 
maps with high velocities (0.75 km s~^ ^ 1.75 km s~^). 
Some infrared H-H objects appear to be associated with 
the B3 component (see Fig. [TT|) . The position- velocity 
diagram indicates that the blueshifted component B8 is 
likely to be the outflow component that appears to move 
away toward south. The blueshifted components B14 
and B7 are located near the edge of Box 1. Although B7 
appears to move from north west to south east from its 
morphology, the driving source remains unknown. There 
is no redshifted counterpart associated with B7. The B7 
component is also associated with some extended green 
objects suggesting that the ambient gas is shocked by the 
B7 component (see Fig. [T^ . The B14 component also 
do not have a redshifted counterpart. There is a com- 
pact dust continuum s o urce immediately south of B14, 
where iBontemps et al.l (12010) found two protostar can- 
didates by the Herschel observations. The Spitzer image 
also shows two compact sources at the same places as the 
dust continuum source and the Herschel protostar can- 
didates (Fig. [T2)) . We suggest that the northern source 
is likely to be the driving source of B14. This interpre- 
tation appears to be supported by the fact that the 1.1 
mm emission is weak at the position of the B14 com- 
ponent, apparently breaking the main fllament into two 
parts (see Fig. [T5k). and bending the northern part of 
the fllament toward south. There are some infrared H-H 
objects in the southern part of B14. These objects may 
originate from B14. 

The relatively-faint blueshifted component B6 is 
clearly seen in the channel maps in the range of 3 km s~^ 
to 5 km s~^, apparently moving toward the south-west 
direction. The possi ble driving source i s one o f the Class I 
sources identifled bv lGutermuth et all (|2008[ ) (source 3 in 
Fig. [TT]) because some H-H objects appear to be ejected 
from it. The B6 component appears to point toward the 
relatively-strong blueshifted component B9. However, it 
is unclear whether B9 originates from the outflow lobe. 
The peak of B9 is about 0.4 pc away from the 1.1 mm 
peak. The position- velocity diagram may imply that this 
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component has different LSR velocity (Vlsr ~ 2 km s~^) 
from the central dense part and does not show the veloc- 
ity structure typical of the outflow lobe (see Fig. [T7b). 
However, there are no YSOs associated with the same 
area as the B9 component. Therefore, this component 
might be a fossil of the outflow lobe or infalling gas to- 
ward the central dense part along a sub- filament. Fur- 
ther investigation is needed to clarify the origin of this 
component. 

There is a redshifted component Rl between the two 
blueshifted components Bl and B2. This redshifted com- 
ponent is clearly seen in the peak intensity map and ap- 
pears to move away from the dense part. This component 
appears to be associated with a blueshifted component, 
which can be seen in the velocity channel maps with 4.75 
km s~^ and 5.75 km s~^. Therefore, this component may 
be almost parallel to the plane-of-sky direction. The pos- 
sible driving source is indicated by a diamond in Fig. Illb 
(source 1), from which many H-H objects are apparently 
ejected. 

The strong redshifted components R2 and R3 appear 
to move away from the central dense part, clearly show- 
ing the velocity structure typical of the outfiow compo- 
nents. Several infrared H-H objects are associated with 
these components. Interestingly, at the area between R2 
and R3, the 1.1 mm emission shows a double-peak fea- 
ture and the 1.1 mm emission tends to be weak on the 
line connecting between R2 and R3 (see Fig. [TSb). In 
addition, there is an area having relatively strong 1.1 mm 
emission just at the north of R2. This morphology leads 
us to the following possibility: the redshifted lobes R2 
and R3 come from the same source and R2 went through 
the dense part, reaching and hitting the dense part lo- 
cated at the north. The R3 component appears to coin- 
cide with the high-velocity component identified by the 
HCO+ (J = 4-3) map (see Fig. [HI). Such high- velocity 
dense gas may suggest that the outfiow is very young. 
As mentioned above, the possible driving source is a 3 
mm source recently discovered by the NMA observations, 
whose position almost coincides with the 1.1 mm peak 
(source 2 in Fig. [11] Fukuda et al. in prep.). However, 
other unknown sources might be the real driving sources 
of these outflow lobes because this central part is very 
dense and crowded. 

There is a faint redshifted component R4 in the south 
west part of the dense part. Several infrared H-H objects 
are associated with the R4 component. This component 
is clearly seen in the channel maps with 11.75 through 
13.75 km s~^, apparently moving toward the south-west 
direction. In fact, the most distant infrared H-H ob- 
ject (K2 in Fig. [T3b ) has a bow shape, consistent with 
the interpretation that it is moving away from the dense 
part. Since a faint blueshifted component labeled by B5 
is associated with the R4 component, this component 
might be almost parallel to the plane-of-sky. The possi- 
ble dTivin£_sourcejs one of the Class I objects identified 
bv lGuterniuth et all (j2008D (source 4 in Fig. [TT]) because 
some infrared H-H objects are apparently ejected from it. 

Several bipolar outfiows can be identified clearly in the 
CO (J — 3—2) data. For example, the southernmost out- 
flow with the B15 and R8 lobes looks relatively powerful 
and clearly shows a bipolar feature. The driving source is 
probably one of the Herschel ca ndidate protostars iden- 
tified bv .Bontemps et al.l (|201Cl[) (see their Fig. 5) whose 



position is (RA, Dec) ~ (18h30m2s, -2°10'20"). The 
candidate protostar is also seen in the 1.1 mm contin- 
uum image presented in Fig. 115b and the Spitzer IRAC 
image in Fig. [TOl as a point-like source. Another bipo- 
lar outfiow can be seen near the upper right corner of 
panel (c) of Fig. [9] The driving source may be either a 
Herschel candidate protostar or a Class object located 
at the upper-right corner of Fig. 5 of iBontemps et all 
()201C1[) . In the northern box with about 13' x 13' area 
(Box 3) we found three molecular outfiows, all of which 
show clear bipolar features: BIO and R5, Bll and R6, 
and B12 and R7. Since this region is out of our Spitzer 
image, we cannot identify the driving sources. 

In Box 2, we identified four additional high- velocity 
components labeled by B13 and R9. These components 
possibly originate from some outfiows. However, further 
observations would be needed to clearly identify the out- 
flow lobes. 

We note that the directions of the outfiow lobes of- 
ten appear to be roughly along the global magnetic 
field direction that i s running almost per pendicularly to 
the main filament (jSugitani et al.ll201ir ). e.g., B2, B5, 
B6, B9, B15, B16, R4, R8, and RIO. However, sev- 
eral components, e.g., Bl, B3, B4, R2, and R3, ap- 
pear to be independent of the global magnetic field di- 
rection. This indicates that the central dense part is 
magnetically supercritical and therefore, the local super- 
sonic turbulence have presumably perturbed the mag- 
netic field significantly. In fact, the recent numeri- 
cal simulations of cluster formation indicate that even 
the moderately-strong magnetic field causes a significant 
misalignment between the global field direction and the 
outflow propa gation directions because of the strong tur- 
bulence Pakamura fc Lill2011| ). 

3.3. Derivation of the Physical Quantities 
3.3.1. Outflow Parameters 

Figure [9] shows that the high- velocity components are 
crowded in the central dense part where the protocluster 
resides and therefore it is very difficult to clearly dis- 
criminate the individual outflow lobes from the data. 
Therefore, we postpone to accurately evaluate the physi- 
cal properties of the individual outfiows until high spatial 
resolution data taken by interferometric observations be- 
come available. In the following, we calculate some global 
physical properties of the outfiowing gas located at the 
cluster center. On the other hand, five bipolar outflows 
mentioned in the previous subsection can be easily iden- 
tified in the CO {J — 3 — 2) map, because they are found 
in relative isolation. However, we postpone to calculate 
the individual outflow parameters of these bipolar out- 
flows because the rms noise level is relatively high with 
a significant scanning effect, precluding accurate estima- 
tion of the outflow parameters. 

Examination of the spectra in the cloud shows that 
ambient ^^CO spectra away from the outfiows appear to 
cover a velocity range of 6 to 11 km s~^. We thus assume 
that the emission in this velocity range include ambient 
emission for the entire cloud. Therefore, the integration 
ranges used are —15 to 4 km s~^ for the blue-shifted 
emission and 11 to 30km s^^ for the red-shifted emis- 
sion, respectively. The outflow masses were calculated 
under the assumption of LTE, a distance of 260 pc, and 
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optically-thin emission in the line wings with an exci- 
tation temperature of 30K. We note that the physical 
quantities estimated below are insensitive to the excita- 
tion temperature. The estimated quantities increase only 
by 20 % ove r the range of T^^ = 20 - 50 K. 

Following iNakamura et al.l ()2011[ ). the outflow mass 
-^^out, the outflow momentum Pout, and kinetic energy 
£^out are estimated as 



and 



-Pout = E^^jl^J - ■^B.V-I/COS^ , 



Eont = E l^^oiVj - V^ys f/ cos' ^ 



(1) 



(2) 



(3) 



respectively, where Mj is the mass of the j-th channel, 
Vj is the LSR velocity of the j-th channel and Vsys is the 
systemic velocity of the driving source. Here, we adopt 
V^ys = 7.5 km s~^. The angle ^ is the inclination angle of 
an outflow, which is generally uncerta i n. He re, we adopt 
^ = 57.3°, following iBontemps et all ()1996[) . The mass 
of the J-th channel, Mj, is expressed as 



i 
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1 fE.TX,.^,^^ 



/ \ arcsec / 
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where i is the grid index on the j-th channel, Nco,i,j is 
the CO column density at the z-th grid on the j-th chan- 
nel, 77 (=0.57) is the main beam efficiency of the ASTE 
telescope, toh2 is the mass of a hydrogen molecule, Xco 
is the fractional abundance of CO relative to H2, and the 
mean atomic weight of the gas fig is set to 1.36. We adopt 
Xco = lO"**. fl is the solid angle of the object and D 
is the distance to the object. From these quantities, we 
calculate the characteristic velocity as Vout = -Pout/-^out- 
Table [2] summarizes the outflow parameters derived 
from the CO ( J = 3 — 2) data. The physical quanti- 
ties presented in Table [2] depend on the adopted velocity 
interval for the integration. For example, if we change 
the velocity intervals as —15 to 6 km s^^ (—15 to 2 km 
s""'^) for the blue-shifted emission and 9 to 30km s"^ 
(13 to 30km s""'^ ) for the red-shifted emission, the total 
mass, momentum, and energy change by -t-65% (—50%), 
20% (-30%), and 5 % (-13 %), respectively. This indi- 
cates that the total momentum and energy do not depend 
strongly on the adopted velocity intervals, although the 
total mass does. This is probably due to the fact that 
the outflow components with lower velocities have lower 
momenta and kinetic energies. As mentioned above, the 
distance to Serpens South remains u ncertain. If we adopt 
D = 410 pc (e.g., iDzib et al.|[20lol ). instead of 260 pc, 
then the mass, momentum, and energy increase further 
by a factor of 2.5. Therefore, a main uncertainty on 
the estimation of the outflow parameters comes from the 



adopted distance to the cloud, rather than the adopted 
velocity intervals. 

3.3.2. Mass of Dense Clump 

The HCO^ (J = 4 — 3) line emission has a critical 
density of 10^ cm"^, and thus traces the dense molecular 
gas. Here, we estimate the mass and radius of the central 
dense clump, using the HCO+ ( J = 4 — 3) data (see Fig. 
W)- The mass of the HCO+ (J = 4 - 3) is evaluated 
from 



M-. 



3.1 X 10" 



V0.57/ 



X 



HCO + 



1.2 X 10- 



/ D 



V260pc 



M y 

arcsec / 



K km s~i 



exp 



43 K 



Mq (5) 



From the above expression, the mass and radius of the 
central clump are estimated to be about 80 Mq and 0.1 
pc, respectively, where we as sumed LTE and th e exci- 
tation temperature of 14 K (jAndre et all l2010[ ). The 
mean density is estimated to be -^ciump/(47ri?^j^j^p/3) ~ 
(3 — 7) X 10^ cm~'^. The estimated mean density is smaller 
than the critical density of the HCO-*- ( J = 4 — 3) line 
by an order of magnitude. This may suggest that the 
central clump contains multiple small subfragments that 
are not spatially-resolved in our observation. In fact, re- 
cent numerical simulations of cluster formation suggest 
that the central dense region of a cluster forming clump 
consists of multiple thin filaments that tend to fragment 
into dens er blobs by the destruction due to the outfiow 
feedback ()Li et al.ll2010D . However, it is difficult to ver- 
ify such a clumpy structure only from our observations. 
Higher spatial resolution observations will be needed to 
reveal the density structure of the dense clump. The 
fractional abundance of HCO"*" relative to H2 is assumed 
to be 1.2 X 10~^, which is adopted from iMaruta et all 
(|2010( ). We note that the fractional abundance of the 
HCO+ remains uncertain toward Serpens South. How- 
ever, from the 1.1 mm data, the 1.1 mm mass in the 
same area as that of the HCO+ clump is evaluated to 
be ~ 10^ Mq, which is not far from the mass estimated 
from the HCO+ emission. Therefore, the adopted HCO+ 
abundance is likely to be reasonable. The FWHM one- 
dimensional velocity width of HCO"'' is estimated to be 
/S.V ~ 2.5 — 3.0 km s^^ near the cluster center. The es- 
timated velocity width corresponds to the Mach number 
of about 5. 

4. DISCUSSION 

According to recent theoretical studies, two main sce- 
narios for cluster formation are proposed. One is the 
rapid, dynamical formation model, in which the large- 
scale flow mainly controls the formation of cluster- 
forming clumps. Star formation is considered to com- 
plete rapidly within a few crossi ng times before initial 
turbulence decays significantly (jHartmann fc BurkertI 
I2007D . In this model, the magnetic field is expected to be 
dynamically weak to promote rapid global gravitational 
collapse. To terminate cluster formation, this model en- 
visions the stellar feedback by the initial star burst that 
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disperses the dense gas away from the parent clumps, 
and eventually the clumps are destroyed. 

The second scenario is the quasi-equilibrium model. 
In this model, star formation process is considered to be 
slow a nd continu e at le ast for several dynamical times 
(jTan et al. 200l iLi fc Nakamu ra WM Matzner 2007. 1. 
Because supersonic turbul ence dissipates qu i ckly in a 
turbu lence-crossing time (jStone et al.l Il998t iMac Lowl 
|1999[) , additional turbulent motions should be injected to 
maintain the supersonic turbulence. In cluster-forming 
clumps, protostellar outflow feedback is considered to be 
a plausible way of driving a significant fraction of super- 
sonic turbulence (iLi fc Nakamural l2006t iMatzneii [20071: 
iNakamura fcLil 120071: ICarroll et al.l I2O1O0 . The proto- 
stellar outflow-driven turbulence can maintain the clump 
close to a quasi- virial equilibrium at least for several 
dynamical times (iMatzned 120071 : INakamura fc Lil 120071 : 
ICarroll et al.l 120101 )^ The magnetic field is expected to 
be dynamically important to imp ede rapid globa l grav - 
itational collapse. Very recently, iSugitani et al.l ()2011[ ) 
found that the global magnetic fields associated with the 
Serpens South filament are well ordered, implying that 
the modcratcly-strong magnetic field played an impor- 
tant role in the formation of the filament. The mag- 
netic structure is likely to be consistent with the quasi- 
equilibrium model. 

To further constrain these cluster formation models, 
it is thus important to clarify the dynamical state of 
the clump and the energy budget of the current out- 
flows in this region, because these two cluster formation 
models predict the very different star formation duration 
timescales. In the following, we attempt to estimate the 
injection rate of the turbulent energy due to the outflows 
and the dynamical state of the dense clump. 

4.1. Turbulence Dissipation and Generation 

Following 'Nakamura et all (|2011| ). we estimate the to- 
tal energy injection rate (mechanical luminosity) due to 
the protostellar outflows and energy dissipation rate of 
supersonic turbulence in the dense clump. From Fig. [51 
the mean length of the outflow lobes can be estimated 
to be 0.1 - 0.2 pc, and the mean velocity along the line- 
of-sight direction is about 7 km s~^ (see Table [2]). From 
these values, the mean dynamical time of the outflows, 
tdyn, can be evaluated as (1 — 2) x lO'' years, assuming 
^ — 57.3°. Thus, the total energy injection rate due to 
the outflows is obtained as 



^out ^ -Eout/^dyn ~ (0.5 — 1)L(^ 



• 



(6) 



On the other hand, the energy dissipation rate of super- 
sonic turbulence is estimated to 



-^turb — C 



1/2MAV^ 
Xd/AV 



(7) 



where C(~ 0.33) is the nondimensional coefficient, AV^ 
is the one-dimensional (ID) FWHM velocity width cor- 
rected for the contribution of the thermal pressure as- 
suming the t e mper ature of 14 K (see Equation [3] of 
iMaruta et al.l ()2010[ )). and Ad is the driving scale of su- 
personic turbulence. From the 1100+ ( J = 4 — 3) data, 
we can derive the energy dissipation rate as 



where we adopted the velocity width and the driving 
scale Xd as AV ^ 2.5 — 3.0 km s~-^ and the clump di- 
ameter of 0.2 — 0.4 pc, respectively. Both the energy 
injection rate Lout and energy dissipation rate Ltmh are 
proportional to the distance to the cloud, D. There- 
fore, the relative importance between Lout and Lturb is 
irrespective of the distance uncertainty. The energy in- 
jection rate due to the outflows is somewhat larger than 
the energy dissipation rate. Therefore, we conclude that 
the outflow feedback can contribute significantly to the 
generation of supersonic turbulence in the clump. 

Since we identified about 5 pairs of the outflows in 
the densest part, the mean outflow momentum for a sin- 
gle outflow may be estimated to be about 2 M© km 
a~^. If we assume the median stellar mass of 0.5 Mq, 
then this gives the outflow momentum per unit stel- 
lar mass of about 4 km s~^, corresponding to the non- 
dimensional parameter / of about 0.04 (see also Maury 
et al. 2009 for the case of NGC2264-C), which gauges 
the strength of a outflow, where the wind velocity of 100 
km s^^ is adopted. This valu e is somewhat smaller tha, n 
the fid ucial values adopted by 'Matzner & McKec' (' 2000D . 
iLL&Nakamura (2006), and Nakamura & Li (2007). The 
estimated value may be somewhat underestimated be- 
cause the outflow lobes are assumed to be optically-thin. 
If the effect of the optical depth is taken into account, 
the value of / is presumably larger than the estimated 
value by a factor of a few. Ev en for a small value of 
/ - 0.1, INakamura fc Lil (|2007^ demonstrated that the 
outflow feedback can supply sufficient momentum to dy- 
namically support a parsec-scale cluster-forming clump 
for several free-fall times, by means of the 3D MHD tur- 
bulent simulations. We note that if the distance to Ser- 
pens South is as large as about 900 pc, th e upper limit 
of the distance to W40 (Rodncv fc Reipu rth 2008,), the 
outflow strength / is estimated to be / 0.4. 

4.2. Dynamical State of Dense Clump 

In the following, we investigate the dynamical state of 
the Serpens South clump, by applying the virial analysis. 
The virial equation for a sphere is given by 



15^ 

2W 



2U + W 



(9) 



where the terms, /, U, and W, denote the moment 
of inertia, internal kinetic energy, and gravitational en- 
ergy, respectively. For simplicity, we neglect the surface 
pressure term. A clump is in virial equilibrium when 
2U + W — 0. The terms U and W are expressed as 



U = 



3MAV^ 
16 In 2 



and 



W 
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R 
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(10) 



(11) 



Lturb - 0.1 - 0.3L, 



O I 



(8) 



respectively, where the values $ and $cr are, respectively, 
the magnetic flux penetrating the clump and the critical 
magnetic flux above which the magnetic fleld can sup- 
port the clump against the self-gravity. Here, we adopt 
$ = for simplicity, which means no magnetic support. 
The value a is a dimensionless parameter of order unity 
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which measures the effe cts of a nonuniform or no nspher- 
ical mass distribution ([Bertoldi fc McKeill992[ ). For a 
uniform sphere and a centrally-condensed sphere with 
p DC r~^, a = 3/5 and 1, respectively. For the HCO+ 
clump, the effects of the nonspherical mass distribution 
appear to be small because the aspect ratios are not so 
far from unity (see also Figure 2 of Bertoldi & McKee 
1992). Here, we adopt a — 1 because the clump seems 
cent r ally- condensed . 

Then, the terms U and W are estimated to be 2U = 
270 - 390 Mo km^ s'^ and W = -350 Mq kuP s^^ 
respectively, which yields the virial ratio, the ratio be- 
tween 2U and W, of 0.8 — 1.1. This value increases 
by a factor of about 1.3, if we adop t the m agnetic field 
strength derived by iSugitani et al.l (|2011[ ). The dense 
clump as a whole is close to quasi-virial equilibrium. On 
the other hand, the outflow kinetic energy of the central 
dense part detected in the CO ( J = 3 — 2) emission is de- 
rived as i^out ^ 24 M© km^ s~^, corresponding to about 
7 % the global gravitational energy, where we adopted 
the value corrected for the inclination angle ^ — 57.3°. 
Therefore, the energy input due to the current outflows 
may not contribute to change the global kinetic energy 
significantly. In other words, the current outflow activity 
observed in Serpens South appears not to be enough to 
disperse the whole cluster-forming clump (see also Maury 
et al. 2009 for the case of NGC2264-C). We note that 
the terms U and W increase with and D^, respec- 
tively, and therefore the virial ratio decreases by a factor 
of 1.6 if we adopt D = 410 pc, instead of 260 pc. Since 
this factor and the factor due to the magnetic support 
almost cancel each other out, the dense clump is likely 
to be close to the virial equilibrium. 

The presence of Class II YSOs in this region suggests 
that star formation has lasted over 10^ years, assuming 
the typical Class II timescale of about 10^ yr. However, 
the low fraction of Class II relative to Class I sources 
implies that the ages of the YSOs in this region may 
be youn ger than those iii othe r nearby cluster-forming 
regions (jCutermuth et al.ll2009( ). The age of this region 
may be 10^~^ yr. This corresponds to a few ~ 20 free-fall 
times of the clump assuming the mean density of 5 x 10^ 
cm where 

According to IGutermuth et al.l (|2008D , about 30 Class I 
and II stars are associated with the area where we de- 
tected the strong HCO+ emission. Assuming the median 
stellar mass of 0.5 Mq, the star formation efficiency in the 
dense part is derived as about 15 — 20 %, as reasonably 
high as those of other nearby cluster-forming clumps. 
This corresponds to about 1 — 5 % of the star formation 
rate per free-fall time. This value is in good ag reement 
with that estimated by 'Krumh olz fc TanI (|2007l ) . Since 
the current outflow activity does not provide enough ki- 
netic energy to destroy the clump, subsequent star for- 
mation is expected to continue. This is consistent with 
the quasi-virial equilibrium model for which cluster for- 
mation continues for a few or more free-fall times. 

5. SUMMARY 



We carried out the CO (J = 3 - 2) and HCO+ (J = 

4 — 3) observations toward the nearest cluster-forming 
IRDC, Serpens South, using the ASTE 10 m telescope. 
The main results are summarized as follows. 

1. We found that many outflow components concen- 
trate in the dense part where the protocluster resides. 
Most of these outflow components appear to move away 
from the dense part. In the northern part, we identi- 
fied only three bipolar outflows, indicating that the star 
formation is less active in the northern part. 

2. Most of the outflow components tend to be anti- 
correlated with the dense gas traced by the 1.1 mm con- 
tinuum emission. Furthermore, the propagation direc- 
tions of t he outflows are across the global magnetic field 
direction (!Sug itani et al.l [20 111) . A couple of strong red- 
shifted outflow components appear to interact with the 
dense gas. 

3. The central dense clump detected by the HCO+ 
( J = 4 — 3) emission has a mass of 80 Mq and a radius 
of about 0.1 pc, with an axis ratio of about 2. The av- 
erage density and velocity width are estimated at about 

5 X 10^ cm"*^ and 2.5 — 3 km s~^, respectively. The dense 
clump has very clumpy structures. Some of the subfrag- 
ments appear to coincide with the positions of th e Class 
candidates identifled bv lBontemps et all ()2010D . 

4. We estimated the global physical quantities of the 
outflows. The total outflow mass, momentum, and en- 
ergy seem smaller than those of the Serpens Cloud Core, 
a nearby typical parsec-scale cluster-forming clump, lo- 
cated about 3° north of Serpens South. However, the 
characteristic outflow speed appears somewhat larger 
than that of the Serpens Cloud Core. This may imply 
that the YSO populations of Serpens South are younger 
than those of the Serpens Cloud Core. 

5. The outflow energy injection rate is likely to be 
somewhat larger than the energy dissipation rate of the 
supersonic turbulence, suggesting that the outflow feed- 
back can significantly contribute to the generation of the 
supersonic turbulence in the dense clump. Assuming the 
median stellar mass of 0.5 Mq, the mean outflow mo- 
mentum per unit stellar mass is estimated to be about 
4 km s~^, under the assumption of optically-thin gas. 
This mean outflow momentum corresponds to the non- 
dimensional outflow parameter of / ~ 0.04, which gauges 
the strength of a outflow. This value of / is somewhat 
small er than the fiducial v alues of / = 0.4 adopted by 
Matzner fc McKed (I2000D . iLi fc Nakan^ (pOOg), and 
iNakamura fc Lil (|2007l ). If we take into account the effect 
of the optical depth, the outflow strength / presumably 
increases by a factor of a few. 

6. The total outflow energy appears signiflcantly 
smaller than the global gravitational energy of the dense 
part where the protocluster is located. In other words, 
it may be difficult to destroy the cluster-forming clump 
by the current outflow activity. This may be inconsistent 
with the dynamical model of cluster formation, for which 
the outflow feedback due to the initial star burst is envi- 
sioned to disperse the dense gas from the cluster-forming 
clump. 
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TABLE 1 

Characteristics of Identified CO ( J = 3 — 2) Outflow Lobes 



Name 




counterpart 




driving source'' 




features 


Bl 




unknown 




source 1 




central part 


B2 




unknown 




unknown 




central part 


B3 




unknown 




unknown 




central part 


B4 




R2, R3 




source 2 




central part, strong 


B5 




unknown 




source 4 




central part, related to R4 


B6 




unknown 




source 3 




central part 


B7 




no 




unknown 




southern part 


B8 




R2, R3 ? 




source 2 ? 




central part, a head of B4 ? 


B9 




unknown 




source 3 ? 




a head of B6 ?, outflow or infall gas 


BIO 




R5 




unknown 




northern part, clear bipolar 


Bll 




R6 




unknown 




northern part, clear bipolar 


B12 




R7 




unknown 




northern part, clear bipolar 


B13 




no 




unknown 




eastern part 


B14 




no 




source 5 




southern part 


B15 




RS 




source 6 




southern part, clear bipolar 


B16 




RIO 




Herschel source 




northern part, clear bipolar 


Rl 




unknown 




source 1 




central part 


R2 




B4 




source 2 




central part, strong 


R3 




B4 




source 2 




central part, strong 


R4 




unknown 




source 4 




central part, related to B5 


R5 




BIO 




unknown 




northern part, clear bipolar 


R6 




Bll 




unknown 




northern part, clear bipolar 


R7 




B12 




unknown 




northern part, clear bipolar 


R8 




B15 




source 6 




southern part, clear bipolar 


R9 




no 




unknown 




northern part. 


RIO 




B16 




Herschel source 




northern part, clear bipolar 


Note. 


— See Sectionl3.2lin 


detail. 










^ Possible counterparts are indicated, "no" means that only one 


component is seen. 






Possible driving sources are 


indicated. Sources 1 through 6 are 


indicated in Figs. Illlandll2l 












TABLE 2 










Global outflow proper,ties derived from CO ( J = 3 — 


2) 




Mass (M 


0) 


Momentum (Mq km s 


Energy (M© 


km2 s-2) Vout (km s"') 


Blue 


0.34 




3.7 




25.9 


10.9 


Red 


0.27 




3.9 




38.7 


14.4 


Total 


0.61 




7.6 




64.6 


12.5 



11 



10 



33 



-1'5«W 



P 



-SHOW 




- 1 



{/J 



a 
■-I 



10" 



I. 



W QflT 
Rigbl Ascension [J2Q0<Dj 



Fig. 1. — Our obser vation area overlaid on the column density map derived from the Herschel observations. The image is taken from the 
left panel of Fig. 1 of lAndre et al.l 12010 ). in which the contrast of the filaments was enhanced by a curvelet transform (see the Appendix 
of Andre et al. 2010). The blue triangles and green stars indicate the prestellar cores and candidate Class protostars identified with 
Herschel, respectively. Our observation area is indicated in cyan boxes that are labeled in numbers. 
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Fig. 2. — CO ( J = 3 — 2) total integrated intensity map toward Serpens South in the velocity range from i'lsr = ~2 km s ^ to +15 km 
s~^. The blow-up of the yellow box is presented in the upper panel of Fig. |6] 
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3Cli4Ds 3O«l20s 




Lj 1 ! 1 ^Uoo 

30tIi4O! lIMlBt 30mfts 18h29ilt«)s 
RA (J20m) 

Fig. 3. — (a) CO (J = 3 — 2) peak intensity map toward the Serpens South protocluster. (b) HCO^ (J = 4 — 3) peak intensity map 
toward the same area presented in pa nel (a), (c) 1.1 m m continuum emission map toward the same area presented in panel (a). The 
detail of the image will be described inlGu termut h et al.l 120 11). For all the panels, the circles and crosses indicate the Class I and Class 
II sources, respectively, identified by IGutermut17'et'^Lrif2008f ) . 
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Fig. 4. — Velocity channel maps of the CO ( J = 3 — 2) emission with the velocity intervals of 1 km s ^. The LSR velocity is indicated 
at the top of each panel. 



16 




-2" (T (T 




-](! U ]« 20 3(1 



2U 10 



to +15 km s 
^. The contours 



Fig. 6. — upper panel: CO (J = 3 — 2) total integrated intensity map in the velocity range from iilsr = ~2 km s^ 
toward the yellow box indicated in Fig. [2] The CO ( J = 3 — 2) integrated intensity is shown in color in units of K km s 
start at 5 K km s~^ at an interval of 5 K km s~^. The boxes indicate the areas where CO { J = 3 — 2) spectra take average (see the lower 
panel of Figure |6]l. lower panel: CO (J = 3 — 2) spectra. Each profile is averaged inside each box indicated in Figure [6^. 
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Fig. 7. — upper panel: HCO+ (J = 4 — 3) total integrated intensity map in the velocity range from i>lsr = 2.25 km s^^ to +15.25 km 
toward the Serpens South protocluster. The HCO"*" ( J = 4 — 3) integrated intensity is shown in color in units of K km s~^. The 
contours start at 2 K km s~^ at an interval of 2 K km s~^. lower panel: HCO"*" ( J = 4 — 3) spectra. Each profile is averaged inside each 
box indicated in Figure [3 
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Fig. 8. — Velocity channel maps of the HCO'*' ( J = 4 — 3) emission with the velocity intervals of 0.2 km s ^. The LSR velocity is 
indicated at the top of each panel. The contour levels go up in 0.15 K km step, starting from 0.2 K km s~^. The HCO"*" ( J = 4 — 3) 
channel maps indicate that the dense clump has clum py structure. The positions of the possible driving sources of the outflows are indicated 
by the white diamonds (sources 1 and 2 in Fig. Illb V Sources 1 and 2 are the possible driving sources of the Rl and R3 (and R2) outflow 
lobes, respectively. The latter coincides with the position of the 1.1 mm peak (see Section 13.21 in detail). 



19 




30m40s 30m20s 30mOs 18h29m40s 

RA (J2000) 



Fig. 9. — (a) Molecular outflow lobes identified from CO ( J = 3 — 2) emission toward the Serpens South protocluster (Box 1). The blue 
contours represent blueshifted ^^CO gas and red contours represent redshifted ^^CO gas. The blue and red contour levels go up in 6 K km 
step, starting from 3 K km s~^. The integration ranges are —9.75 to 3.75 km s ~^ for blueshifted gas and 11.25 to 29.25 km s~^ for 
redshifted gas. (b) Same as panel (a) but for the northern part of the Serpens South filament (Box 3). The blue and red contour levels go 
up in 6 K km step, starting from 6 K km . (c) Same as panel (a) but for the southern part of the Serpens South filament (Box 2). 
The blue and red contour levels go up in 6 K km s~^ step, starting from 6 K km s~^. 
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18:^:40 . 18r30J;0^ ' -: - 18:30:00 

■ RiifitAscsnsipn (J2000) 




Fig. 10. — (a) Three-color image of Serpens South with IRAC band 1 (3.6 /im, blue), IRAC band 2 (4.5 ^mm green), and IRAC band 4 
(8.0 fi m, red), (b) Same as panel (a) but the positions of Class I {magenta circle) and II {cyan circle) sources identified bv .Gutcrmuth et al.l 
lj2008l ) are overlaid. The positions of the infrared H-H objects are indicated by the crosses. 




Fig. 11. — (a) Blow-up of the Spitzer IRAC i mage toward the central region of Serpens South. The diamonds indicate the positions of 
the possible outflow driving sources (see Section |3 21 in detail), (b) Same as panel (a) but the blueshifted and redshifted CO (J = 3 — 2) 
out flow lobes are overlaid on the image by the blue and red contours. The contours and circles are the same as those of Fig. |9ji and Fig. 
1101 respectively. 
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Fig. 12. — (a) Blow-up of the Spitzer IRAC image toward the southern part of Serpens South. The s quar es indicate the positions of the 
same objects shown in Fig. 115b . The blow-ups of the areas labeled by Kl and K2 are presented in Fig. 1131 (b) Same as panel (a) but the 
blueshifted and redshifted C O (J = 3 — 2) outflow lobes are overlaid on the image by the blue and red contours. The circles and crosses 
are the same as those of Fig. [10] 
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(a) 





Fig. 13. — Blow-ups of the areas indicated by Kl and K2 in Figs. Illb and 112b ,. An infrared H-H object having a bow shape is seen near 
the center of each panel. 




Fig. 14.— HCO+ (J = 4 - 3) and CO (J = 3 - 2) line profiles toward source 2 indicated in Fig. [TTJ3. The HCO+ { J = 4 - 3) shows a 
clear blue-skewed profile, indicating that the envelope gas is infalling toward the central object. 
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Fig. 15. — (a) Blow-up of the CO ( J = 3 — 2) outflow lobes located at the southern part of Box 2. The gray colors indicate the 1.1 
mm continuum image. The upper and lower diamonds are the possible driving sources of Rl and R3 (and R2), respectively. The latter 
coincides with the position of the 1.1 mm peak, (b) Same as panel (a) but for the northern part of Box 1. The possible driving sources of 
the outflows are indicated by the green circles whose positions are determined from the positions of the sources seen in the Spitzer IRAC 
image. In panels (a) and (b), the blue contours, red contours, and labels are the same as those of Fig. |9] 
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RA (J2000) 

Fig. 16. — High velocity components identified from HCO+ (J = 4 — 3) emission toward the Serpens South protocluster. The color image 
indicates the CO (J = 3 — 2) integrated intensity map in the range of 2 km to 15 km s~^. The white and black contours represent 
blueshifted and redshifted gas detected by HCO+ ( J = 4 — 3), respectively. The white and black contour levels go up in 0.5 K km s~^ and 
1.0 K km s~^ steps, starting from 1.0 K km and 1.0 K km s~l, respectively. The integration ranges are 3.25 to 5.75 km s ~^ fo r th e 
blueshifted gas and 9.75 to 15.25 km for the redshifted gas. The squares indicate the positions of the same objects shown in Fig. 115b . 



